Entanglement induced by nonadiabatic chaos 
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We investigate entanglement between electronic and nuclear degrees of freedom for a model nona- 
diabatic system. We find that entanglement (measured by the von Neumann entropy of the subsys- 
tem for the eigenstates) is large in a statistical sense when the system shows "nonadiabatic chaos" 
behavior which was found in our previous work [Phys. Rev. E 63, 066221 (2001)]. We also discuss 
non-statistical behavior of the eigenstates for the regular cases. 
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Quantum information processing (QIP) is one of hot 
topics in many branches of science 1]. One important 
point is how to implement a quantum computer in real 
systems, and many possibilities have been theoretically 
suggested and experimentally tested. One candidate can 
be molecular systems because highly excited molecules 
have dense quantum states, which can be manipulated by 
laser fields. Some quantum logic gates in such a system 
can be built by using optimal control theory |2| , and are 
actually realized in a molecular system 

In highly excited molecules or laser-driven molecular 
systems, nonadiabatic transition (NT) is a rule rather 
than an exception ^ |^, i.e., we have to consider not 
only electronic or nuclear degrees of freedom (DoF) , but 
both at the same time. In such a case, a fundamental 
issue related to QIP is how much (quantum) entangle- 
ment is produced between electronic and nuclear DoF in 
molecular systems, because entanglement is a key ingre- 
dient in QIP. 

Here we investigate a two-mode- two-state (TMTS) sys- 
tem which has two electronic DoF and two nuclear (vi- 
brational) DoF with a nonadiabatic coupling |3, Q • This 
is a (minimum) model NT system which shows "quan- 
tum chaos" behavior, i.e., statistical properties of energy 
levels and eigenstates are similar to those of a random 
matrix system Q. If many electronic DoF are involved, 
the similar system has a naive classical limit, and its dy- 
namical property of entanglement has been already ad- 
dressed in [3. In such cl CclSG, Oil quantum chaological 
view is effective, and we can say much about a quantum 
system by studying its classical limit 0- However, the 
situation is different and more difficult here, because the 
TMTS system does not have a naive classical limit due 
to discreteness of the electronic DoF and deserves 
further attentions in view of entanglement. 

The TMTS system in the diabatic representation is 
described by the following Hamiltonian: 
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where Tkin is the kinetic energy, Vi {i = A, B) is the 



potential energy for state i defined by 
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X cos 6 — {y — a) sin 0, 
X sin 6 + {y ~ a) cos 9, 
X cos 9 + {y + a) sin 9, 
—X sin 9 + {y + a) cos ( 
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The geometrical meaning of the parameters is shown in 
Fig.n Note that we have just used a harmonic potential 
for each state. 
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FIG. 1: A shematic representation of the TMTS system. The 
distance between the minima of the potential is 2a, and the 
angle between the relevant crossing seam (dotted line) and 
the primary axis of each potential (dashed line) is 8. Inset: 
The perspective of the TMTS system. The potential minima 
are different with Ae = eb — = 0.173. 

Here the Duschinsky angle 9 Tl'l and the nonadiabatic 
coupling constant J are two important parameters for the 
system; The latter induces entanglement between elec- 
tronic and vibrational DoF. We solve this Hamiltonian 
according to the procedure in 6], and obtain the cigen- 
energies and eigen- vectors. The fc-th eigenvector can be 
written as 
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where = {\A} + \B})/V2, |2)) = + \B))/V2, 

\A), \B) are the electronic bases for diabatic surfaces A 
and B, respectively, \j) represents 2D harmonic eigen- 

functions, and |</.f ^) = = Ej (^ff^lj)- 

From this eigenvector, we can construct a reduced den- 
sity operator for the electronic DoF as 
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where c\^^ are actually all real numbers. 

The measure of entanglement we choose here is the von 
Neumann entropy of the subsystem defined by 



= -TrlpWlogpf'^)} 
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where Af ^ (« = 1,2) is an eigenvalue for the 2x2 matrix, 
Eq. ©. A note is in order: the value of the entropy is 
the same if we use the reduced density operator for the 
vibrational DoF. We took the electronic DoF because the 
2x2 matrix is very easy to diagonalize, and to interpret 
the result as shown below. 

First we show the J dependence of the results fixing 
9 = Tr/6: As we can see in Fig. 12 the entropies for the 
case of J = 1.5 assemble around its maximum — 
log 2, whereas those of the other cases (J — 0.3, 7.5) are 
rather broadly distributed. This condition of entangle- 
ment is very similar to that of quantum chaos behavior 
found in [Q]: When both J and 6 have "intermediate" 
values (J ~ 1 and 9 ~ 7r/4), the system shows quan- 
tum chaos behavior, i.e., the nearest neighbor spacing 
distribution becomes the Wigner type, A3 statistics a 
log curve, and the amplitude distribution of the eigen- 
states Gaussian. To further confirm this, we show the 9 
dependence of the results fixing J = 1.5 in Fig.O This 
result also nicely corresponds to the previous condition 
for the quantum chaos behavior. From these results, we 
can conclude that, in a statistical sense, the condition 
for quantum chaos behavior to arise in the TMTS sys- 
tem is very similar to that for entanglement production in 
eigenstates to arise in the same system. This conclusion 
result supports the previous study which uses entan- 
glement production as an indication of quantum chaos 
behavior. (It is also noted that the calculation of entan- 
glement production is rather easier than those like the 
nearest neighbor spacing distibution or A3 statistics be- 
cause there is no fitting procedure.) 

However, as noticed in Fig.|21 the amount of entangle- 
ment strongly varies depending on each cigcnstate for the 
cases of J = 0.3 and 7.5. To consider this problem, we 
rewrite the entropy (entanglement) using the vibrational 



bases for each electronic state: |(/>^^) and {'^''b')- These 
states are connected to the above states \(t>[''^) and {(f)'^'') 




Energy level number: k 

FIG. 2; J dependence of entanglement production measured 
by the von Neumann entropy as a function of the energy level 
number. Top: J = 7.5. Middle: J = 1.5. Bottom: J = 0.3. 
The Duschinsky angle is fixed as 6 = ir/G. 
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Energy level number: k 

FIG. 3: 9 dependence of entanglement production measured 
by the von Neumann entropy as a function of the energy level 
number. Top: 6 = tv/3. Middle: 6 = tv/6. Bottom: 6 = 0.0. 
The nonadiabatic coupling is fixed as J = 1.5. 
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FIG. 4: S^^' (+), lAP^'^l (*), and 15^*^1 (x) as a function of 
energy level number, (a) J = 0.3 (weakly nonadiabatic case), 
(b) J — 7.5 (strongly nonadiabatic case). (Lines are just for 
guiding eyes.) The Duschinsky angle is fixed as 9 = n /6. 



Hence the density operator, Eq. lO, is represented as 
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where we have introduced two new parameters: Sj^g is 
the overlap between the fc-th eigenstates on surface A 
and B, and AP^^ is the half of the population difference 
between the fc-th eigenstates on surface A and B. Us- 
ing these parameters, the eigenvalues for the entropy is 
written as 

x[''}^\±^/\s^ln^\A^^. (16) 

From this relation, for the entropy to be large, both |>5'^'^^| 
and |AP^*^| should be small. We can numerically confirm 
this for the strongly "chaotic" case {J = 1.5), and this 
property might be derived from the random matrix the- 
ory. For less "chaotic" cases (J — 0.3,7.5), the situation 
is different: As shown in Fig. ^ there is a strong corre- 
lation between S^'^ and jAP^'^l for the weakly nonadia- 
batic case (J = 0.3), whereas between S^'^ and \s'^^\ for 
the strongly nonadiabatic case (J — 7.5). On the other 
hand, is'xll ^ for the former and \AP^^\ ~ for the 
latter. This is interpreted as follows: For the former, the 
eigenstates "reside" on diabatic surfaces A and B which 
are tilted each other. Thus the overlapping between the 

eigenstates \s'^g\ becomes small because the nodal pat- 
terns for the eigenstates are also tilted (see Fig. 5 in |0|). 
For the latter, the eigenstates "reside" on adiabatic sur- 
faces, and the amplitudes of them on diabatic surfaces A 
and B are similar (see Fig. 7 in [^), hence jAP^^I ~ 0. 

Let us focus on the weakly nonadiabatic case (J = 0.3). 
In the range of fc = 800 to 820, the lowest entangled state 
is 805-th, and the highest is 815-th [Fig. H (a)] . In Figs. El 
andEl we show the two eigenstates on diabatic surfaces A 
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FIG. 5: A less entangled regular state: 805-th eigenstates on 
diabatic surfaces A (A) and B (B) for the weakly nonadiabatic 
case: J = 0.3. Pi{x,y) = \{x,y\(l)i)f {i = A,B). Note that 
the scale for (B) is smaller than that for (A). The Duschinsky 
angle is fixed as ^ = 7r/6. 
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FIG. 6: A strongly entangled regular state: 815-th eigen- 
states on diabatic surfaces A (A) and B (B) for the weakly 
nonadiabatic case: J = 0.3. Pi{x,y) = \{x,y\(f>i}f {i — A,B). 
The Duschinsky angle is fixed as 6 = n /6. 



and B. As anticipated from the above argument, there is 
a large population difference on surface A and B for the 
less entangled state, whereas there is not for the strongly 
entangled state. The latter situation means that even a 
regular state can strongly entangle. Note that, albeit we 
desymmetrized the system with a finite Ae = — es, 
we have this entangled state for the regular case. (If we 
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do not desymmetrize the system, i.e., Ae = 0, we easily 
have entangled states for both regular and chaotic cases 
because of the symmetry.) Thus we must be cautious to 
use the entanglement production as a manifestation of 
quantum chaotic behavior 9J, though it is fine to use it 
in a statistical sense. 

In this paper, we investigated quantum entanglement 
production between electronic and nuclear (vibrational) 
degrees of freedom for a nonadiabatic system. We found 



that the condition of the entanglement in the eigenstates 
to appear is statistically very similar to that of quantum 
chaos behavior to show up. We also discussed the non- 
statistical behavior of the entanglement production and 
interpreted it using the eigenstate properties on diabatic 
surfaces. It will be interesting to analyze other nonadi- 
abatic systems like Jahn- Teller molecules IT^ in light of 
entanglement production. 



[1] M.A. Nielsen and I.L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

[2] J. P. Palao and R. Kosloff, Phys. Rev. Lett. 89, 188301 

(2002) ; Phys. Rev. A 68, 062308 (2003). 

[3] J. Vala, Z. Amitay, B. Zhang, S.R. Leone, and R. KoslofF, 
Phys. Rev. A 66, 062316 (2002). 

[4] H. Nakamura, Nonadiabatic Transition: Concepts, Ba- 
sic Theories and Applications (World Scientific, Singa- 
pore, 2002); C. Zhu, Y. Teranishi, and H. Nakamura, 
Adv. Chem. Phys. 117, 127 (2001). 

[5] K. Takatsuka, Y. Arasaki, K. Wang, and V. McKoy, Fara- 
day Discuss. 115, 1 (2000); Y. Arasaki, K. Takatsuka, 
K. Wang, and V. McKoy, Phys. Rev. Lett. 90, 248303 

(2003) . 

[6] H. Fujisaki and K. Takatsuka, Phys. Rev. E 63, 
066221 (2001). See also H. Fujisaki and K. Takatsuka, 
J. Chem. Phys. 114, 3497 (2001); H. Higuchi and 
K. Takatsuka, Phys. Rev. E 66, 035203(R) (2002). 



[7] E.J. Heller, J. Chem. Phys. 92, 1718 (1990). 
[8] M. C. Gutzwiller, Chaos in Classical and Quantum Me- 
chanics (Springer- Verlag, New York, 1990). 
[9] K. Furuya, M.C. Nemes, and G.Q. Pellegrino, 
Phys. Rev. Lett. 80, 5524 (1998). 

[10] Note, on the other hand, that there is a sound correspon- 
dence between the TMTS system and its mapping (quasi- 
classical) system. See H. Fujisaki, |quant-ph / 0401 136] 
Phys. Rev. E (in press). 

[11] See J. Tang, M.T. Lee, and S.H. Lin, J. Chem. Phys. 119, 
7188 (2003), and references therein. 

[12] H. Koppel, W. Domcke, and L.S. Cederbaum, 
Adv. Chem. Phys. 57, 59 (1984); D.M. Leitner, 
H. Koppel, and L.S. Cederbaum, J. Chem. Phys. 104, 
434 (1996); H. Yamasaki, Y. Natsume, A. Terai, 
and K. Nakamura, Phys. Rev. E 68, 046201 (2003); 
A.P. Hines, CM. Dawson, R.H. McKenzie, and G.J. Mil- 
burn, quant-ph/0402016J 



